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Catalyst-driven electrolysis of water is considered as a "cleanest" way for hydrogen production. Finding 
cheap and abundant catalysts is critical to the large-scale implementation of the technology. 
Two-dimensional metal dichalcogenides nanostructures have attracted increasing attention because of their 
catalytic performances in water electrolysis. In this work, we systematically investigate the hydrogen 
evolution reduction of metal dichalcogenides monolayers based on density- functional-theory calculations. 
We find that metal disulfide monolayers show better catalytic performance on hydrogen production than 
other metal dichalcogenides. We show that their hydrogen evolution reduction strongly depends on the 
hydrogen coverage and the catalytic performance reduces with the increment of coverage because of 
hydrogenation-induced lower conductivity. We further show that the catalytic performance of vanadium 
disulfide monolayer is comparable to that of Pt at lower hydrogen coverage and the performance at higher 
coverage can be improved by hybridizing with conducting nanomaterials to enhance conductivity. These 
metal disulfide monolayers with lower overpotentials may apply to water electrolysis for hydrogen 
production. 



Energy has been a key in the development of every sphere of human society. The traditional usable energy 
sources, including fossil fuels and coal, will fall short of the demand of sustainable development over the 
long term, and their continued use produces harmful side effects such as pollution that threatens human 
health and greenhouse gases associated with climate change, which have triggered considerable world-wide effort 
to explore renewable green energy alternatives 1 . Hydrogen is considered to be one of the most important 
candidates because of its abundance and clean and renewable nature. As an ideal clean energy carrier for future, 
hydrogen can be produced from a variety of energy resources, has the highest energy density per unit mass, and 
produce the least polluting since it can be extracted from natural resource such as water or biomass and its use 
produces water only. The electrolysis of water is considered as a well-known principle to produce oxygen and 
hydrogen gas in a sustainable fashion 2 " 6 . The key component in electrochemical reduction of water is the catalyst 
for hydrogen evolution reduction (HER). The well-known catalysts in the electrolysis of water are noble metals, 
such as platinum, due to their superior electrocatalytic properties 7 " 10 . The application of noble metal catalysts in 
large-scale production of hydrogen is limited by their high cost and low abundance, although lots of efforts had 
been done by tuning the composition of the catalyst to modify their electronic structures 11 " 24 . On the other hands, 
considerable efforts have been carried out to search alternative catalysts with lower cost and abundance 22 " 25 . 
Basically, an advanced catalyst for the enhanced electrochemical hydrogen evolution reaction should reduce the 
HER reaction overpotential and consequently increase the HER efficiency 3 " 6 ' 25 . 

The transition-metal dichalcogenides with the formula of MX 2 , where M is a transition metal element from 
group IV (Ti, Zr, or Hf), group V (for instance V, Nb or Ta) or group VI (Mo or W), and X is a chalcogen (S, Se or 
Te), have attracted increasing attention for their applications in electrolysis of water 26 " 43 . These materials have 
crystal structures consisting of weakly coupled sandwich layers X-M-X, where one M-atom layer is enclosed 
within two X layers and the atoms in layers are hexagonally packed 43 . The overall symmetry of transition metal 
dichalcogenides can be hexagonal or rhombohedral, and the metal atoms have octahedral (IT) or trigonal 
prismatic (2H) coordination 43 . Experimental and theoretical studies had suggested that the electrocatalytical 
activity of 2H-MoS 2 in electrolysis of water is contributed to its edges 27 " 29 , which are metallic if they are zigzag 44,45 . 
Recently, Voiry et al. 37 reported that metallic 1T-WS 2 nanosheets showed better HER performance than semi- 
conducting 2H-WS 2 , which can be further improved by strain engineering. As compared with the noble metals, 
we may be sure that only MX 2 monolayers with high conductivity, and MX 2 nanoribbons/nanoparticles with 
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Figure 1 | The representative structures of pristine MX 2 (a) and fully H-covered MX 2 (b) top and (c) side views. 



metallic edges can promise the excellent HER activity. To date, a 
comprehensive study on the HER performances of MX 2 monolayers 
as well as its origin has not been available. In this work, the applica- 
tions of MX 2 monolayers in electrolysis of water are systematically 
investigated based on the calculations of density-functional theory 
(DFT). We predict that the HER performances of MX 2 monolayers 
depend on M, X and H- coverage. We show that VS 2 is comparable to 
Pt for electrolysis of water at lower H- coverage and its catalytic 
activity can be further enhanced by improving conductivity. 

Results and discussion 

In our calculations, we focus on metallic/semimetallic 2H transition- 
metal dichalcogenides, MX 2 (M = Nb, Ta, and V; X = S, Se, and Te) 
because 2H phase is more stable and high conductivity is essential to 
the electrolysis of water. The MX 2 unit cells with trigonal prismatic 
(2H) coordination are first optimized to obtain the lattice para- 
meters. The optimized structures of MX 2 (M = Nb, and Ta; X = 
S, Se, and Te) (see Supporting Data, Table I) from our calculations 
are consistent with the reported experimental and theoretical data 46 . 
We see that the lattice parameters of NbX 2 are almost equivalent to 
those of TaX 2 , while are larger than those of VX 2 (Supporting Data, 
Table I). To investigate the HER performances of MX 2 monolayers at 
various hydrogen coverage on their surfaces, the geometries of MX 2 
monolayers with one surface fully covered by hydrogen atoms (MX 2 - 
H) are relaxed to find out the effects of hydrogen coverage on their 
lattice parameters. The hydrogen atoms are adsorbed on the top of X 
atoms (Figure 1), where is the most stable position 37,4748 . The opti- 
mized structures of MX 2 -H (Table II in Supporting Data & Figure 2) 
show that the lattice constants (a) are extended by 2.0 to 3.8%, 1.7 to 
3.3%, and 3.4 to 4.5% for NbX 2 , TaX 2 , and VX 2 , respectively, where 
the extension increases as X changes for S to Se, further to Te 
(Figure 2). The thicknesses (c) and the X-M bonds of the monolayers 
are reduced by the H-coverage on their surfaces (Figure 2). The 
calculated bond lengths are about 1.4, 1.5 and 1.7 A for S-H, Se-H, 
and Te-H bonding, respectively (Supporting Data, Table II). 

According to the Sabatier principle, the optimal catalytic activity 
of material for HER can be achieved on a catalytic surface with 
intermediate binding energies (or free energies of adsorption) for 
reactive intermediates 15 , which can be quantified by analyzing the 
reaction free energy of hydrogen adsorption (AG H ) 13,15 ' 37 ' 49 . The 
optimum value is around AG H = 0. To obtain the reaction free 
energy, we calculate AG H for various H coverage on MX 2 monolayers 
as following: 

AG H = AE H + AE ZPE - TAS H ( 1 ) 



where AE H is the hydrogen chemisorption energy defined as: 
AE H = E(MX 2 + nH) - E(MX 2 + {n-l)H)- ^E(H 2 ) 



(2) 



where n is the number of H atoms adsorbed on a MX 2 monolayer. 
Full coverage refers to one hydrogen atom per X atom adsorbed on 
one side of the MX 2 monolayer. The AG H as a function of the hydro- 
gen coverage can be obtained by changing n. E(MX 2 + nH), E(MX 2 ) 
and E(H 2 ) in Eq. (2) are the energies of monolayer with hydrogen 
atoms (n), pure MX 2 monolayer (n = 1), and hydrogen molecule, 
respectively. AS H is the difference in entropy. The entropy of adsorp- 



tion of 1/2 H 2 is AS H = — 1/2S° H2 , where S° H2 is the entropy of H 2 in 
the gas phase at standard conditions. AE ZPE is the difference in zero 
point energy between the adsorbed and the gas phase, related to the 
reaction: l/2H 2 (g) — > H* 9 where H* denotes a hydrogen atom 
adsorbed on the surface. AE Z pe ~ TAS H is about 0.24 e V 13 ' 37,49 , sim- 
plifying Eq. (1) to AG H = AE H + 0.24. 

Rectangle supercell (Figure 1) is employed to study the different 

i 

hydrogen coverage on MX 2 monolayer, -(z = l~6). The full cov- 

6 

erage (- or 1) is that all of X atoms on one side of the monolayer 
6 

adsorb hydrogen atoms (Figures lb&lc). A supercell with double 
size in y direction of Figure 1 is used for — coverage, that is, one 

H atom is attached to one of the 12 X atoms. — coverage is that one H 

6 

atom is attached to one of the 6 X atoms (Figure 1). For other 

coverage (-(/ = 2~5)), all of the possible H-adsorption configura- 
6 

tions on the six X atoms are considered, where the adsorption energy 
is calculated by two methods - "average" of all configurations and 
"most stable" configuration. We first consider the supercell (p- 
supercell) based on monolayer pristine lattice parameters (from 
Table I in Supporting Data). The calculated adsorption energies 
show that the overpotentials are positive, and AG H increases with 
the increase of H coverage on the MX 2 monolayer (Figure 3), which is 
in contrast with that of Pt (AG H trends on zero with the increase of H 
coverage) 1315 ' 49 . We see that the calculated AG H follows the same 
trend regardless of the calculation methods (Figures 3a, 3c, and 3e 
for "average" method, and Figures 3b, 3d, and 3f for "most stable" 
method) and the difference on values induced by the method is 
minor. Importantly, the A G H of MS 2 is lower than those of MSe 2 
and MTe 2 by 40% at the same H coverage (Figure 3), indicating that 
the HER performances of MS 2 monolayers are better than those of 

MSe 2 and MTe 2 monolayers. The overpotentials at — coverage from 

"average" method (Figures 3a, 3c, & 3e) are 0.061, 0.192, and 

0.007 eV for NbS 2 , TaS 2 , and VS 2 , respectively. At ^ coverage, the 

AG H for NbS 2 , TaS 2 , and VS 2 are 0.1 12, 0.266, and 0.051 eV, respect- 
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Figure 2 | The difference of lattice parameters, including lattice constant 
(a), monolayer's thickness (c), and X-M bond length (X-M), between 
pristine and H-covered MX 2 monolayers. 
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Figure 3 | The calculated overpotentials as a function of H-coverage of MX 2 monolayers in p-supercell by "average" method: (a), (c), and (e); and by 
"most stable" method: (b), (d) and (f). 



ively. We see that the HER performances of NbS 2 , TaS 2 , and VS 2 
monolayers are better than that of 1T-WS 2 monolayer because their 
AG H at relatively high H-coverage is lower than that of 1T-WS 2 at 

lower H-coverage (AG H = 0.28 eV at — coverage and AG H = 

16 

2 

0.36 eV at — coverage) 37 . The HER performance of VS 2 at a cov- 
16 

erage up to - is also better than that of MoS 2 edges (AG H = 
6 

0.08 eV) 28 . For comparison, we compose the curves of AG H of MS 2 
(M = Nb, Ta, and V) monolayers as a function of H-coverage into 
Figure 4a. We see that the HER performances of MS 2 monolayers 
increase as M changes from Ta to Nb, further to V. The overpotential 

of VS 2 at — coverage is 10 and 30 times lower than those of NbS 2 and 
12 6 

TaS 2 . The VS 2 -based catalyst for electrolysis of water may be com- 
parable to that of Pt at lower H-coverage (up to — coverage) because 

of its near-zero overpotential (0.007 eV). The difference at HER 
performance, however, decreases with the increase of H-coverage 

on the monolayer. At full coverage (- or 1), VS 2 only take advantage 

6 

of NbS 2 and TaS 2 by 6 and 20%, respectively. 

We also calculate the hydrogen adsorption energy in the supercell 
(h-supercell) constructed on the lattice parameters of fully H-cov- 
ered monolayer (from Table II in Supporting Data). We see that the 
HER trend as a function of H-coverage (Figure 5) is the same as that 
in p-supercell. The MS 2 monolayers still show the best HER perform- 
ance at same H-coverage because their overpotentials are smaller 
than those of MSe 2 and MTe 2 monolayers (Figure 5). For all of 
MX 2 monolayers, AG H in h-supercell is reduced because of the 
extended lattice constants (that is equivalent to the effect of strain 
in Ref. 37). The difference of AG H between p-supercell and h- 



supercell, however, becomes smaller as the H-coverage increases, 
indicating that the strain shows less effect on HER of the monolayer 

if the H-coverage is high. At — coverage, the overpotentials from 

"average" method (Figures 5a, 5c, & 5e) are —0.07, 0.098, and 

—0.159 eV for NbS 2 , TaS 2 , and VS 2 , respectively, and at - coverage, 

6 

the corresponding AG H are 0.004, 0.174, and —0.138 eV, respect- 
ively, which should imply that NbS 2 shows the best HER perform- 
ance at lover H-coverage in h-supercell. With the increase of 
H-coverage on the MS 2 monolayers' surfaces, AG H increases 
(Figure 4b). After 25% H-coverage, VS 2 gives out the best HER 
performance because its overpotential is positive and lower than 
those of NbS 2 and TaS 2 (Figure 4b). By carefully comparing the 
energies of systems calculated from p-supercell and h-supercell, we 
find that the effect of H-coverage on lattice parameters is negligible 
2 

when it is less than the 33% (- coverage). Therefore, we combine the 
6 

calculated AG H of MS 2 from p-supercell for H-coverage up to - and 

6 

2 6 

h-supercell for H-coverage from - to - into Figure 6. Clearly, we see 

6 6 

that VS 2 monolayer shows the best HER performance in all of con- 
sidered systems. 

To investigate the possible origin of HER performance and the 
effects of H-coverage, the partial densities of states (PDOSs) of M, X, 
and H are calculated. The calculated PDOSs clearly show the change 
of conductivities of MS 2 monolayers with the H-coverage (Figure 7 
and supporting data SI ~ S5). We see that the pristine VS 2 mono- 
layer is metallic, where the Fermi level is within the middle energy 
band and near to its bottom (Figure 7a). The middle energy states are 
dominated by V d electrons. With introducing H atoms to the mono- 
layer's surface, the energy sates shift down to lower energy region, or 
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Figure 4 | The calculated overpotentials as a function of H-coverage of MS 2 monolayers by "average" method in p-supercell (a) and in h-supercell (b). 



the Fermi level shifts up to high energy region within the middle 
band (Figure 7b). The Fermi level further shifts up with increasing H- 
3 

coverage (Figure 7c ~ h). At - coverage, the Fermi level of the system 
6 

is close to the middle energy band top. Corresponding to semi- 
conductor's band structure, we may state the middle band as valence 
band and the band above as conduction band. The up-shift of Fermi 
level to the middle band (valence band) top leads to the reduction of 
conductivity. This is further illustrated on VS 2 monolayer with full 
H-coverage, where the system is an intrinsic semiconductor or insu- 
lator because its Fermi level is within the band gap (Figure 7h). The 
PDOSs of NbS 2 and TaS 2 in p-supercell and h-supercell, and VS 2 in 
h-supercell reveal the same origin of the reduction of performance 

1.6 



under high H-coverage (Supporting data, SI ~ S5). The analysis of 
PDOSs shows that the conductivities of the MX 2 monolayers are 
reduced by the H-coverage, especially after half H-coverage on the 
surfaces, leading to the decreases of their HER performances 
(Figures 3 — 6). This finding may also apply to the MX 2 nano- 
structures with metallic edges because their conductivities should 
be reduced when edges are fully saturated by hydrogen atoms 44 . 
This is different from Pt-related catalysts, whose conductivities have 
not been affected by the H-coverage on their surfaces 13,15 ' 49 . How to 
improve the conductivities of MS 2 -based catalysts under high H- 
coverage should be a challenge for their application in electrolysis 
of water. One of ways is to mix these MS 2 monolayers with high 
conducting nanostructures, such as graphene 33,50 . 
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Figure 5 | The calculated overpotentials as a function of H-coverage of MX 2 monolayers in h-supercell by "average" method: (a), (c), and (e); and by 
"most stable" method: (b), (d) and (f). 
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Figure 6 | The combined overpotentials as a function of H-coverage of 

2 

MS 2 monolayers by "average" method in p-supercell (H-coverage < - ) 

2 6 
and in h-supercell (H-coverage > -). 

6 

For comparison, the HER performance of MoS 2 monolayer is 
calculated. We see that its overpotential (AG H ) (Supporting data, 

56) is larger than those of MX 2 (M = Nb, Ta, and V; X = S, Se, 
and Te) monolayers considered in this work (Figure 6). At lower and 

full H-coverage (— ,- , and 1), its overpotential is about 2.0 — 2.1 eV. 
12 6 

At moderate H-coverage, its overpotential is larger than 1.1 eV. The 
calculated densities of states show that pure and H- covered MoS 2 
monolayers are intrinsic or n-type semiconductors (Supporting data, 

57) . These results confirm that the surface of MoS 2 monolayer is inert 
to electrolysis of water, and the most active sites are at the edges or 
defects of MoS 2 nanostructures 27 29 . However, our study show that 
the surfaces of MX 2 (M = Nb, Ta, and V; X = S, Se, and Te) mono- 
layers, especially VS 2 , are very active, which further enhance their 
ability for electrolysis of water due to higher contacting area with 
water. 

Conclusions 

The DFT-based first-principles calculations are carried out to invest- 
igate the hydrogen evolution reduction of MX 2 monolayers. We find 
that MS 2 monolayers are better than MSe 2 and MTe 2 monolayers in 
electrolysis of water, especially VS 2 , which shows the best HER per- 
formance because of its lower overpotential. We show that the HER 
performances of MX 2 monolayers strongly depend on the H- 
coverage on their surfaces. With increasing H-coverage, the perform- 
ance is reduced because of the reduction of conductivity. We also 
show that the strain may improve the HER performance at relatively 
low H-coverage. We further predict that their HER applications at 
high H-coverage can also be achieved by improving their conduct- 
ivities, such as hybridization with metallic nanostructures. It is 
expected that the MS 2 monolayers, especially VS 2 , may find applica- 
tions to electrolysis of water for hydrogen production. 

Methods 

The first-principles calculations are carried out to investigate the hydrogen evolution 
reduction of transition-metal dichalcogenide monolayers. Our calculations are based 
on the density functional theory (DFT) 51 and the Perdew-Burke-Eznerhof generalized 
gradient approximation (PBE-GGA) 52 . The projector augmented wave (PAW) 
scheme 53 ' 54 as incorporated in the Vienna ab initio simulation package (VASP) 55 is 
used in the study. The Monkhorst and Pack scheme of k point sampling is used for 
integration over the first Brillouin zone 56 . A15X 15X 1 grid for k-point sampling for 
geometry optimization of unit cells, and an energy cut-off of 450 eV are consistently 
used in our calculations. A sufficiently large supercell is used so that the monolayers in 
neighbouring cells in the vertical direction are separated by a vacuum region of at least 
20 A. Good convergence is obtained with these parameters and the total energy was 
converged to 2.0 X 10" 5 eV/atom. The error bar (or uncertainty) of the DFT 
calculation is less than 5 meV. The thermodynamic processes via Tafel pathway are 




-3-2-10123 



Energy (eV) 

Figure 7 | The calculated partial densities of states of pristine (a) and 
various H-covered (b ~ h) VS 2 monolayers in p-supercell. 

calculated 17 . The effect of solvent on the HER performance of VS 2 monolayer is 
investigated by including H 2 0 molecules in the systems with various H- 
coverage 18 ' 21 . 

The calculated thermodynamc processes of HER on VS 2 monolayer via Tafel 
channels and the effect of solvent on overpotentials are included in Supporting Data. 
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